Phase slips in superconducting films with constrictions 
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A system of two coplanar superconducting films seamlessly connected by a bridge is studied. We 
observe two distinct resistive transitions as the temperature is reduced. The first one, occurring 
in the films, shows some properties of the Berezinskii-Kosterlitz-Thouless (BKT) transition. The 
second apparent transition (which is in fact a crossover) is related to freezing out of thermally acti- 
vated phase slips (TAPS) localized on the bridge. We also propose a powerful indirect experimental 
method allowing an extraction of the sample's zero-bias resistance from high-current-bias measure- 
ments. Using direct and indirect measurements, we determined the resistance -R(r) of the bridges 
within a range of eleven orders of magnitude. Over such broad range, the resistance follows a simple 
relation R{T) = Rn exp [-{c/t){l - t)^/'^], where c = AF{0)/kTc is the normalized free energy 
of a phase slip at zero temperature, t = T/Tc is normalized temperature, and Rn is the normal 
resistance of the bridge. 



I. INTRODUCTION 

Thermally activated vortex-like excitations (topologi- 
cal defects) of the superconducting condensate is the pri- 
mary source of dissipation in mesoscopic superconduct- 
ing structures^ These fluctuations take different forms 
in one-dimensional (ID) and two-dimensional (2D) sys- 
tems. In 2D thin films the fluctuations are known 
to be broken vortex-antivortex pairsSiSiliSiSiLSiSiiSiii 
while in ID wires the resistance is due to phase 
Slipsii^i"ii^ii5ii6,i7,i8,i9,20,2i important difference 

between these two types of fluctuations is that vortices 
and antivortices form bound pairs below a certain crit- 
ical temperature, known as the Berezinskii-Kosterlitz- 
Thouless (BKT) transition temperature, while phase 
slips and anti-phase-slips are unbound at any finite tem- 
perature. Thus the resistance of ID wires is greater than 
zero at any finite temperature due to the presence of 
phase slips, which are described by the theory of Langer, 
Ambegaokar, McCumber, and Halperin (LAMH)!i2ii^ 

Here we report a study of structures in which both 
types of fluctuations can coexist, namely thin films con- 
taining constrictions, which are comparable in size to the 
coherence length. The goals of this work are (i) to test 
the applicability of the LAMH theory for short and rather 
wide constrictions and (ii) to test the effect of the vortex- 
antivortex sea existing in the thin film banks adjacent to 
the constriction on the phase slippage rate on the con- 
striction itself. For this purpose we fabricate and measure 
a series of thin superconducting MoGe films22i (which 
are about 15 fim wide) interrupted by constrictions or 
"bridges" (see Fig. ^). The width of the narrowest 
point of the bridges is in the range of 13-28 nm, i.e. a 
few times larger than the coherence length (we estimate 
^(0) w 7 nm for our MoGe films^'^). Two resistive transi- 
tions are observed in such samples indicating that vortex- 
antivortex pair binding-unbinding transition (if any) and 
thermally activated phase slip processes occur separately. 
For T > Tbkt the contribution of vortex-antivortex pairs 
is dominant. On the other hand, below Tbkt the trans- 
port properties are determined by the phase slip pro- 



cess on the bridge, which may be regarded as a vortex- 
antivortex pair breaking assisted by the bridge. 

Using direct and indirect techniques we have tracked 
the sample's resistance within a range of eleven orders 
of magnitude. The resulting R{T) curves are compared 
with the LAMH theory. Regardless of the large width 
of the bridges and their shortness, the shape of the mea- 
sured R{T) curves is in perfect agreement with the overall 
shape of the curves computed using the standard LAMH 
theory (note that this theory was originally derived for 
very long wires that are much thinner than the coher- 
ence length). The only disagreement found with LAMH 
is that the pre-exponential factor had to be modified in 
order to obtain a reasonably low critical temperature of 
the bridges. (The critical temperature is used as an ad- 
justable parameter in the fitting procedure.) Following 
the argument of Littlsi^ we arrive at the conclusion that 
the pre-exponential factor should be simply Rn and ob- 
tain a good agreement with measured curves. The mea- 
surements show that the bridges with intermediate di- 
mensions (i) allow phase slippage which does not quench 
at any finite temperature, (ii) behave independently of 
the thin film banks, and (iii) exhibit a higher rate of 
phase slippage in the cases when the width of the bridge 
is smaller and therefore when the coupling between the 
thin film banks is weaker. 

Before presenting our experimental results we give a 
brief summary of the BKT theory of topological phase 
transitions and the LAMH theory of thermally activated 
phase slips. In thin superconducting films, even in the 
absence of a magnetic field, an equal population of free 
vortices and antivortices is expected to occur. The BKT 
theory predicts a universal jump in the film superfiuid 
density at the characteristic temperature Tbkt, lower 
than the mean field critical temperature of the film Tcq. 
Such a jump is related to the vortex-antivortex pair bind- 
ing through a logarithmic interaction potential between 
free vortices^^^. Applied currents can break bound pairs 
producing free vortices and leading to non-linear V{I) 
curves. Above Tbkt the linear resistance of a film is 
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given by the Halperin-Nelson (HN) formula&S 



-2v/6(rco-rBKT)/(r-r] 



BKTJ 



^HN = lO.SbRnj exp 

(1) 

where i?„.t is the normal state resistance per square of the 
film and 5 is a non-universal constant. Note that the HN 
equation predicts zero resistance for temperatures below 
the BKT phase transition temperature Tbkt- 

The LAMH thcorjii*i24iii4 applies to narrow supercon- 
ducting channels, in which thermal fluctuations can cause 
phase slips, i.e. jumps by 27r of the phase difference of 
the superconducting order parameter. In unbiased sam- 
ples the number of phase slips (which change the phase 
difference by 2tt) equals the number of anti-phase-slips 
(which change the phase difference by —2tt). An applied 
bias current pushes the system away from the equilib- 
rium and the number of phase slips becomes larger than 
the number of anti-phase-slips. Thus a net voltage ap- 
pears on the sample, which can be calculated, following 
LAMH, as = S0/2e (below we will also discuss an al- 
ternative approach to the voltage definition). Here h is 
Planck's constant, e is the electron charge, and (p is the 
rate of change of the phase difference between the ends 
of the wire. During the phase slip process the energy of 
the system increases since the order parameter becomes 
suppressed to zero in the center of the phase slip. Ther- 
mal activations of the system over this free energy bar- 
rier AF(T) occur at a rate given by (17(T)/27r)e-'^-^/'=^. 
If the bias current is not zero, then the net rate of the 
phase slippage is </> = f^(r)(e-'^^+(^)/'=^^ - e^'^^-'^^y'''^). 
Here / is the bias current, and AF^ and AF^ are the 
barriers for phase slips and anti-phase-slips correspond- 
ingly (these two barriers become equal to each other at 
zero bias current). The attempt frequency derived from 
a time-dependent Ginzburg-Landau (GL) theory, for the 
case of a long and thin wire, isi^ 



n(r) 



L 1 



/AF 



1/2 



(2) 



where T is the temperature of the wire, and L/£^{T) is the 
length of the wire measured in units of the GL coherence 
length ^(T). The attempt frequency is inversely propor- 
tional to the relaxation time tql = Trh/8k{Tc — T) of the 
time-dependent GL theory, with Tc being the mean field 
critical temperature of the wire (or of the bridge, as in 

1/2 

our discussions below). The factor (Ai^/fcT) ' provides 
a correction for the overlap of fluctuations at different 
places of the wire and the factor L/£^{T) gives the num- 
ber of statistically independent regions in the wire^^. The 
free energy barrier for a single phase slip is give n^?i^? by 



8%/2i?2(T) 



(3) 



which is essentially the condensation energy den- 
sity H^{T)/8tt multiplied by the effective volume 



8\/2A^(r)/3 of a phase slip {A is the cross-section area 
of the wire). 

A bias current / causes a non-zero voltage (time aver- 
aged) given by 



V 



^ffie-^^/'=^sinh(///o) 



(4) 



where Iq = 'iekT/h {Iq = 13.3 nA at T = 1 K). Differen- 
tiation of this expression with respect to the bias current 
/ gives the differential resistance 



dV/dl 



e/n 



cosh(///o) 



(5) 



The dependence of the attempt frequency and free energy 
on the bias current is neglected in this derivation. In the 
limit of low currents I <ti Iq, Ohm's law is recovered 



el. 



kT J 



(6) 

where Rq = h/{2e)'^ = 6.5 kfl. In this approach the fluc- 
tuation resistance does not have any explicit dependence 
on the normal resistance of the wire. 



II. EXPERIMENTAL SETUP 

The sample geometry is shown schematically in Fig. 
The fabrication is performed starting with a Si wafer 
covered with Si02 and SiN films. A suspended SiN bridge 
is formed using electron beam lithography, reactive ion 
etching, and HF wet etchingiSi The bridge and the en- 
tire substrate are then sputter-coated with amorphous 
Mo79Ge2i superconducting alloy, topped with a 2 nm 
overlayer of Si for protections. The resulting bridges 
are 100 nm long with a minimum width w « 13-28 nm 
as measured with a scanning electron microscope (SEM) 
(Fig. njj). All samples are listed in Table I. 



80 nm 




SiN supporting 
membrane 



FIG. 1: (a) Sample schematic. MoGe film (black) of thick- 
ness d = 2.5-4.5 nm is deposited over a SiN membrane (gray) 
substrate with a constriction of width w. (b) An SEM mi- 
crograph of a typical sample. The MoGe coated SiN bridge 
(gray) is suspended over a deep trench (black). 



Transport measurements are performed in a pumped 
"'He cryostat equipped with a set of rf-filtered leads. The 
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linear resistance R{T) is determined from the low-bias 
slope (the bias current is in the range of 1-10 nA) of the 
voltage versus current curves. The high-bias differential 
resistance is measured using an ac excitation on top of a 
dc current offset generated by a low-distortion function 
generator (SRS-DS360) connected in series with a 1 Mil 
resistor. One sample was measured down to the mfi level 
using a low temperature transformer manufactured by 
Cambridge Magnetic Refrigeration. 



III. RESULTS 

First we compare a sample with a hyperbolic constric- 
tion ("bridge sample") with a reference sample, which 
is a plain MoGe film of the same thickness, without 
any constriction ("film sample"). Both are fabricated on 
the same substrate simultaneously. A resistive transition 
measured on the film sample is shown in Fig. [51 The 
HN fit generated by Eq. ^ is shown as a solid line and 
exhibits a good agreement with the data, yielding a BKT 
transition temperature of Tbkt — 4.8 K and the mean 
field critical temperature Tco = 4.91 K. Such good fit sug- 
gests that the transition observed in the banks might be 
the BKT transition, although a more extensive set of ex- 
periments is necessary in order to prove this assumption 
rigorously. As expected, Tbkt is slightly lower than Tcq. 
The inset of Fig. |2] compares the i?(T) measurements 
of the "film" (open circles) and the "bridge" (solid Hne) 
samples. At T = 4.8 K the R{T) curve for the film sam- 
ple crosses the i?, = axis with a nonzero (and large) 
slope, in agreement with the behavior predicted by the 
HN resistance equation (Q. Nevertheless, unlike the film 
sample, the bridge sample shows a non-zero resistance 
even below the BKT transition temperature predicted 
by Eq. Such resistive tails, occurring at T < Tbkt, 
have been found in all samples with constrictions. 

In Fig. Othe R{T) curves for five samples with bridges 
are plotted in a log- linear format. The resistance of sam- 
ple Bl has been measured down to the mQ range using 
a low temperature transformer. Two resistive transitions 
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Sample 


w (nm) 


Rn (n) 


Tc (K) 


no (K) 


d (nm) 


P 


Al 


27 ±4 


1380 


3.88 


3.90 


2.5 


1.47 


Bl 


13 ±4 


1650 


4.80 


4.91 


3.5 


0.723 


B2 


28 ±4 


1320 


4.81 


4.91 


3.5 


1.10 


CI 


13 ±4 


1440 


5.16 


5.50 


4.5 


0.653 


C2 


27 ±4 


680 


5.39 


5.50 


4.5 


2.21 



TABLE L Sample parameters, including the width of the 
constriction (w), determined from SEM images, normal resis- 
tance of the bridge (Rn), determined from the R{T) curves 
(at a temperature slightly below the resistive transition of the 
thin film banks), critical temperature (Tc), determined from 
^WL (T) fits given by Eq. ©) , critical temperature of the film 
(Tco), film thickness (d), and a geometrical fitting parameter 
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FIG. 2: Low-bias resistance versus temperature dependence 
(open circles), measured on a thin film (d = 3.5 nm) without 
constriction. The solid line is a fit to the Halperin-Nelson 
theory (Eq. Q). (Inset): Resistance of the film without 
constriction (multiplied by a constant factor), shown as open 
circles, is compared to the sample with a hyperbolic bridge 
(B2), shown by the solid line. The only qualitative difference 
is the presence of a "resistive tail", observed on all samples 
with constrictions. 



are seen in each curve as the temperature decreases. The 
first transition is the superconducting transition in the 
thin film banks adjacent to the bridge. The second tran- 
sition corresponds to the resistive tail mentioned above. 
In order to understand the origin of the second transition 
it should be compared to the LAMH theory. 



IV. DISCUSSION 

Below we analyze the resistive tails found on samples 
with constrictions and demonstrate that they are caused 
by the phase slip events localized on the bridges and be- 
have independently of the adjacent thin film banks. The 
analysis indicates that no BKT (no vortex-antivortex 
binding within the constrictions) or any other type of 
transition occurs on the constrictions and that the phase 
slips and anti-phase-slips are unpaired at any nonzero 
temperature due to thermal fluctuations. This is demon- 
strated below by fitting the R{T) curves with the LAMH- 
like fitting curves. 



A. LAMH attempt frequency for a short bridge 

In order to compare our results to the LAMH we have 
to take into account the small length of the bridge, which 
does not allow more than one phase slip at a give time. 
Therefore the attempt frequency fl{T) of Eq. ^ can 
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be simplified. First, it has a term L/£_{T) that accounts 
for the number of independent sites where a phase sUp 
can occuni^. Since each of our samples has only one 
narrow region where phase slip events can happen, we 
take L/^(T) = 1. Second, the coefhcient {AF/kTy/^ 
which takes into account possible overlaps of phase slips 
at different places along the wire^ is taken to be unity 
also. This is because for short hyperbolic bridges (not 
much longer than the coherence length) it is reasonable 
to expect that there is only one spot, i.e. the narrowest 
point of the bridge, where phase slips occur. As a result, 
we obtain the attempt frequency for a short hyperbolic 
bridge flwL = I/tgl (the abbreviation "WL" stands for 
"weak link"). This attempt frequency can be combined 
with the usual form of the LAMH resistance in Eq. © 
and can be used to fit the experimental R(T) curves (be- 
low the resistive transition of the films). Although such 
fits follow the data very well, there is one inconsistency 
that is they require the critical temperature of the bridge 
to be chosen higher than the critical temperature of the 
films, which is unphysical for such system. We attempt 
to modify the pre-exponential factor in order to resolve 
this inconsistency, as discussed below. 



a 

DC 




T(K) 

FIG. 3; Low-bias resistance for five different samples with 
bridges. The parameters of the samples are given in Tabled] 
The data points are shown by open symbols. Solid lines are 
fits to the "bridge" phase slip model given by Eqs. ||7J and 
0. 



B. Modification of the prefactor 

Since the exact expression is unknown, we approximate 
the resistance of a constriction (weak link) as 



The exponential factor here is that of the LAMH the- 
ory and the prefactor is simply the normal resistance 
of the bridge. This expression (Eq. {T))) can be justi- 
fied by the following argument: the duration of a single 
phase slip (i.e. the time it takes for the order parame- 
ter to recover) is ^ tql and the number of phase slips 
occurring per second is ~ ilwi,{T) exp [— AF(T)/fcr], 
with the attempt frequency being the inverse GL relax- 
ation time f^wL — ^/tgl, as was argued above. There- 
fore the time fraction during which the constriction is 
experiencing a phase slip (i.e. when superconductiv- 
ity is suppressed on the bridge) is the product of these 
two values, i.e. / = (rGL)(l/TGL) exp[-AF{T)/kT] = 
exp [—AF{T)/kT]. Following Littleji^ it can be assumed 
that the bridge has the normal resistance Rn during the 
time when a phase slip is present (i.e. when the bridge 
is in the normal state), and the resistance is zero other- 
wise (when there is no phase slip). Thus we arrive at the 
averaged resistance for a bridge or a small size weak link 
R^^ ^ f X Rn + (1 - f) X = Rn exp [~AF/kT] as 
in Eq. ^ . Note that unlike in the LAMH theory, in the 
present formulation the fiuctuation resistance is directly 
linked to to the normal state resistance of the sample. 

In order to compare Eq. Q to the experimental re- 
sults, an explicit expression for the energy barrier A^Vl 
for a phase slip localized on the bridge is required. Start- 
ing with the usual form.^^ derived for a long ID wire and 
some well known results from BCS and GL theoryjiii^ 
we find that AFwl(O) = (8\/2/3)(iJ2(0)/87r)A^(0) = 
0.83kTcRqL/ Rn£,{0) where L is the length of the wire. 
Using Rn — PnL/A, the free energy barrier for a weak 
link is 



AFwl(O) = OMkTc 



PwdRq 

Pnm 



(8) 



where w is the width of the bridge, d is the film thickness, 
Pn is the normal resistivity, and A = wd. The parameter 
(3 measures the ratio of the phase slip length along the 
bridge to the effective length of a phase slip in a ID 
wire, which is equal to 8\/2^(T)/3. Finally, assuming 
the same temperature dependence of the barrier as in 
the LAMH theory, i.e. AF{T) = AF(0)(1 - T/Tcf^^, 
we arrive at the expression for the bridge fluctuation- 
induced resistance: 



^WL (T) = Rn exp 



-0.83 



PwdRq 

Pnm 



3/2 rp 

1-^1 4 



T 



Rwl{T) = Rn e-^^^-/'^. 



(7) 



(9) 

The fits generated by Eq. are shown in Fig. |2| 
as solid lines. An impressively good agreement is found 
for all five samples. In particular, sample Bl measured 
using the low-temperature transformer, shows an agree- 
ment with the predicted resistance i?wL over about seven 
orders of magnitude, down to a temperature that is more 
than two times lower compared to the critical tempera- 
ture of the sample. Only two fitting parameters are used: 
(3 and Tc (listed in Table QJ. The other parameters re- 
quired in Eq. (Q, including R^, d, w, ^(0), and p„ « 180 
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0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 



I(HA) 

FIG. 4: Differential resistance as a function of the dc bias 
current for sample B2. Experimental data are denoted 
by open symbols and the solid lines are fits to dV/dl — 
ii(T) cosh(///o). Temperatures from left to right are 4.12, 
3.92, 3.80, 3.64, 3.45, 3.36, 3.26, 3.16, 3.07, 2.80 and 2.68 K. 

liQ. cm are known liSiiSiS The fits give quite reasonable 
values for the critical temperature of the bridges, in the 
sense that they are slightly lower than the corresponding 
critical temperatures of thin films of the same thickness, 
as expected. This fact supports the vahdity of Eq. 0. 
Such good agreement also indicates that the dissipation 
in a thin film with a constriction at T < Tbkt is solely 
due to thermal activation of phase slips on constrictions. 
As expected, /3 « 1 for all samples and the larger f3 values 
are found on wider constrictions. 



C. Determination of the linear resistance from 
high bias current measurements 

We now discuss the non-linear properties of films with 
constrictions. Measurements of the differential resistance 
versus bias current, dV/dl vs. /, are plotted in Fig. 0]on 
log-linear scale. Using these results it is possible to dis- 
tinguish between the BKT mechanism, which leads to a 
power-law V{I) dependence, and the phase slippage pro- 
cess, which is characterized by an exponential V{I) de- 
pendence (Eqs. Q and ©). From Fig. 2|it is clear that 
at T < Tbkt and sufficiently low currents the dependence 
of the differential resistance on bias current is exponential 
(it appears linear on the log- linear plots). Thus it is ap- 
propriate to compare the results with the LAMH theory. 
Equation Q can be written as dV/dl = R{T) cosh(///o), 
where R(T) is the temperature-dependent zero-bias re- 
sistance. Using this relation, we fit the differential re- 
sistance data and use R{T) as a fitting parameter, as 
shown in Fig. ^ by solid lines, each corresponding to a 




T(K) 

FIG. 5: Resistance vs. temperature curve for sample B2. 
Open circles represent data that have been directly mea- 
sured while filled boxes give the resistance values determined 
by fitting the dV/dl curves of Fig. |1] using the formula 
dV/dl — i?(r) cosh(7//o). The solid and the dashed curves 
give the best fits generated by the Rwh{T) {Tc = 4.81 K) and 
i?LAMH(r) (Tc = 5.38 K) formulas, respectively. 



fixed temperaturei^i The fitting procedure illustrated in 
Fig. 21 gives us a powerful indirect method of determina- 
tion of the zero-bias resistance (it is implicitly assumed 
that the ratio of the rates of thermally activated and 
quantum phase slips (if any) is independent of the bias 
current). This method is useful when the temperature is 
low and the resistance of the sample is below the reso- 
lution limit of the experimental setup. Thus, by fitting 
the dV{I) I dl curves, we obtained the zero-bias resistance 
R{T) down to very low values 10~^ fl). This method 
was systematically applied on sample B2 and the results 
are shown in Fig. [31 as solid squares. The open circles 
in Fig. |31 represent the zero-bias resistance obtained by 
direct measurements at low bias currents. The two sets 
of data are consistent with each other. The solid curve in 
Fig. I^lis a i?wL fit obtained using Eq. ©. An excellent 
agreement is seen in a wide range of resistances spanning 
eleven orders of magnitude. This re-confirms that the 
thermally activated phase slip mechanism is dominant in 
the bridge samplesSS for T < Tbkt- We emphasize that 
the critical temperature of the bridge, which is used as 
an adjustable parameter, is found to be Tc = 4.81 K. As 
expected, the Tc of the bridge is slightly lower than the 
critical temperature of the film electrodes Tco = 4.91 K. 

The usual LAMH expression i?LAMH (Eq. ©), which 
applies to thin superconducting wires jiSiiSiSSiS^ can also 
be used to fit our data. The overall shape of the fitting 
curve (dashed curve in Fig. [SJ agrees with the data as 
well as with the i?wL fit- The drawback of the usual 
LAMH formula is that the critical temperature of the 
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bridge, which is used as an adjustable parameter, turns 
out considerably higher than the film transition temper- 
ature. For example, the dashed line fit in Fig. 13 is gen- 
erated using Tc — 5.38 K which is larger than the film 
critical temperature T^q — 4.91 K. This apparent en- 
hancement of the critical temperature of the bridge must 
be an artifact, because a reduction of the dimensions of 
MoGe samples always leads to a reduction of the critical 
temperaturei^ On the other hand, the Tc extracted from 
the fits made using Eq. @ are almost equal and slightly 
lower than the film T^o (Table I), as expected. 

A rapid decrease of the LAMH resistance at tempera- 
tures very close to the critical temperature refiects the be- 
havior of the LAMH attempt frequency which approaches 
zero as T ^ Tc. The LAMH resistance is proportional to 
the attempt frequency so we observe i? — > as T — s- Tc 
(dashed curve in Fig. ) . Such behavior is unphysical 
and occurs since the LAMH theory is not applicable very 
near T^. It should be emphasized that some of our mea- 
sured bridges are wider than ^(0), yet the thermally acti- 
vated phase slip model agrees well with the data. This is 
in agreement with the prediction (Ref. 0, p. 510) that 
superconducting channels of width w < 4.4^(T) should 
exhibit a ID behavior, i.e. nucleation of vortices is un- 
favorable in such channels. Such condition is true for all 
of our samples. 

V. SUMMARY 

Fluctuation effects in thin films interrupted by "hy- 
perbolic" constrictions is studied. The measurements 



show two separate resistive transitions. The higher- 
temperature transition shows some properties of a BKT 
transition in the films (follows the HN formulae). The 
second apparent resistive transition is explained by a 
continuous reduction of the rate of thermally activated 
phase slips with decreasing temperature. A quantitative 
description of the fluctuation resistance of narrow and 
short superconducting constrictions is achieved. For this 
purpose we have modify the LAMH expression for the 
resistance of a one-dimensional nanowire. An indirect 
method that enables us to trace the resistance variation 
over eleven orders of magnitude is suggested, based on 
the analysis of the nonlinear effects occurring at high bias 
currents. The phase slippage model is found applicable 
in the entire range of measured resistances, suggesting 
that quantum phase slipst^^ do not occur in this samples, 
in the studied temperature interval, which extends below 
Tc/2 for one sample (Bl)). 
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